


f 


j! S nN 
iby 4 ‘ ts Me , 









Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1971 


Temperature distribution in thin plates welded underwate 


Staub, James Andrew. 


Massachusetts Institute of Technology 


http://ndl.handle.net/10945/24871 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
_ (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist | et Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
TT \ KNOX appointed -- and published — scholarly author. 
http://www.nps.edu/library 






LIBRARY 


Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 


TEMPERATURE DISTRIBUTION IN THIN PLATES 
WELDED UNDERWATER 


James Andrews Staub 











c00S-fr 


TEMPERATURE DISTRIBUTION IN THIN PLATES 
WELDED UNDERWATER 
by 


JAMES ANDREW STAUB, JR. 
1 
/} 


A.B. MATHEMATICS, ST. VINCENT COLLEGE 
(1963) 


Beoell.5., PENNSYLVANIA STATE UNIVERSITY 
(1964) 


SUBMITTED IN PARTIAL FULFILLMENT OF THE 
REQUIREMENTS FOR THE DEGREE OF 
NAVAL ENGINEER AND THE DEGREE OF 
MASTER OF SCIENCE IN MECHANICAL ENGINEERING 
at the 


MASSACHUSETTS INSTITUTE OF TECHNOLOGY 


SUNE LOT EL 


SIGNATURE OF AUTI 


CERELPUEDUBY i 


BCCECTED ay <a... 


ACTUAL LS PL St, TER Sh VA LN Ned ee as ek te as 


Be Nilpnediity: 
TOOHOS auvA 


AYVUEIT XONN AIG 


— a a ae ee ee 





LIBRARY 
NAVAL P&STGRADUATE SCHOOB 
MONTEREY, CALIF. 93940 


TEMPERATURE DISTRIBUTION IN THIN PLATES 


WELDED UNDERWATER 


By 


JAMES ANDREW STAUB, JR. 


Submitted in partial fulfillment of the requirements for 
the degrees of Naval Engineer and Master of Science in Mechani- 
cal Engineering at the Massachusetts Institute of Technology, 
June, 1971. 


ABSTRACT 


During underwater (wet) welding the metal within the heat- 
affected-zone experiences a rapid fluctuation in temperature 
which results ina loss of ductility and a susceptibility to 
stress-corrosion cracking. The object of this study was to 
develop and verify a model which would predict the temperature 
distribution such that the metallurgical and mechanical pro- 
perties of the underwater weld could be determined. 


The eguation of heat flow in two-dimensions with source 
terms was converted to its finite difference equivalent and 
solved by succesSive iterations with a digital computer. The 
solution incorporated all non-linear terms in the original 
equation and accounted for the variation of the thermal pro- 
perties of the metal with temperature. The interior boundary 
in the temperature field was found to have a definite shape 
whth was a function of the important welding parameters. The 
Shape function was correlated by experiment using methods . 
established for normal surface welding. Heat losses were as- 
sumed to occur by boiling and radiation heat transfer. 


The predicted and measured temperature distributions were 
compared. It was found that the model predicted peak tempera- 
tures with reasonable accuracy but that comparison of cooling 
curves showed considerable disagreement. Since the predicted 
cooling curve presented a sharper slope than the measured value 
it was considered that the assumption of conventional boiling 
heat transfer processes was invalid. 


Thesis Advisor: Prof. K. Masubuchi 
Title: Associete Professor of Naval Architecture 
and Marine Engineering 
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CHAPTER I 


INTRODUCTION 


The object of this thesis was to develop a model of an 
underwater welding process such that the temperature history 
of the metal could be predicted. Particular attention was 
given to temperature changes within the heat affected zone. 
The model was constructed with consideration for a) the varia- 
tion with temperature of the thermal properties of the base 
metal, b) a finite heat source of predictable configuration, 
c) the heat losses to the liquid environment, and da) the non- 
linear components of the heat flow equation. Three major 
areas within the literature were investigated: a) under- 
water, b) normal surface welding, and c) boiling heat trans- 


fer. 


This study was conducted in recognition of the poten- 
tial benefits that could be realized if underwater welding 
were a viable process. Predicted trends in ocean exploration 
and exploitation indicate the use of large structures that 
could only be fabricated or repaired on site within the 
liquid environment. Even today, considerable savings could 
be realized if ship hulls, off-shore drilling platforms, buoys 
and salvage work were constructed or repaired without ex- 
plicit dependence on the requirement that the welding opera- 
tion be conducted in air. Unfortunately, underwater welds 


are not structurally sound due to reduced ductility, and 
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susceptibility to stress corrosion and hydrogen cracking. 


A reasonably accurate model, one that predicts peak 
temperatures and cooling rates, was seen as a necessity to 
examine the metallurgical changes occurring in the weld 


zone as function of welding parameter. 


A. UNDERWATER WELDING: STATE-OF-THE-ART 


Underwater welding, as a oer is best defined by 
considering its two subcategories: wet and dry welding. 
Dry underwater welding is a welding procedure in which the 
liguid environment is excluded from the vicinity of the work- 
piece by a structure such as a caisson. The structure is 
large enough to contain one or more workers plus the neces- 
sary welding equipment. The atmosphere within the chamber is 
usually at depth pressure and may, or may not, consist of 
life-supporting gases. This system's utility is restricted 
to symmetrical workpieces about which the chamber can be 
fitted and sealed. Hence, dry underwater welding is most 
used by oil companies laying undersea pipeline or repairing 
the submerged supporting structure on off-shore drilling 
rigs. Essentially, this procedure differs from normal air 
welding only by the pressure of the atmosphere surrounding 


the work. 


Wet welding occurs when the liguid environment is al- 


lowed to interact directly with the joining process: the 


ae 
ae eee | — a 


CE 
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most notable result BSing Bi esyonlal quenching cf the liquid 
metal. The utility of this process is not ae by 

the geometry of the workpiece, but rather by the quality of 

the product. No wet weld is considered safe when used in 

load bearing structures. The eee quench produces a weld that 
has 50% of the ductility of air welds, and is subject to 
extreme stress corrosion cracking and hydrogen embrittle- 

ment. The object of the thesis, then, was to construct, 

and verify, a model that might predict the temperature cy- 


cles in a metal as it is being wet welded. 


The necessity for such a nedmieis twofold: a) considera- 
ble time and expense could be saved if current ocean struc- 
tures could be repaired or:.fabricated in situ, and b) future 
projections of man's involvement with the ocean require the 
use of structures that can only be fabricated, by reason of 
their size, within the ocean environment. Currently, the 
most economical means of joining metal is by welding, and it 
would be extremely expensive if each new process, each new 
structural material, with their multitude of variables, had 
to be investigated individually by actual experimentation. 
Unfortunately, the state of wet welding has generally stag- 


nated since its inception during WW II. 


From its development in the 1940's until today, the 
most widely used wet welding process is shielded metal-arc. 


This is the same method used in air welding and is commonly 
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femgerredsto as “covered electrode" or "stick electrode” 
welding. Generally, wet welding has merely borrowed what- 
ever equipment was already fewela nec eee air welding. The 
only adaptation to the new environment being a coat of water- 
proof material for the electrode and an electrically-insulated 
holder. Only recently have the British begun to develop 
electrodes specifically designed for use under water. How- 
ever, the details of the program are considered proprietary 


information. 


The basic research conducted thus far has been an ex- 
amination of the mechanical snamemeies of the flux and the 
completed weld. Thus, Waugh and Eberlein’? investigated 
PitsecOVverings tO Cetermine Ehose most Suitable for under— 
water application.* Work at the Naval Engineering Station~° 
produced some factors of safety that must be applied when 
wet welding. Hibshman and Jensen® evaluated the properties 
of welds in low-carbon steel, under water and in air. The 
data revealed that ductility was markedly reduced. As far 
as it is known, however, no attempt has been made, either on 


a theoretical or experimental basis, to determine the tem- 


perature history of the metal during the wet welding process. 


B. NORMAL SURFACE WELDING 


A review of the literature for air welding, however, 


indicated that considerable effort has been devoted to 


* Numbers refer to listing in Bibliography. 





- 
predicting the temperature cycles during the welding process. 
In 193% pememenai-° obtained an exact solution to a modi- 
fied heat flow equation. He siemens his model by neglec- 
ting non-linear terms, treating the conductivity as constant, 
and assuming a point source. However, the resulting equation 
for temperature distribution contained all the important weld- 
ing parameters: weld speed, heat input, plate thickness, and 
the thermal properties of the metal. Furthermore, the equa- 
tion was in quasi-steady state form so that the location of 
the point source was at the origin of a stationary coordinate 
system. Experimental verification by Rosenthal?® and 
Schmorber indicated that the model was reasonably accurate 

at some distance from the source of heat, but in adjacent 
areas (within the heat affected zone, or haz) the degree of 
error was unacceptably large. The primary cause for the 
discrepancy was due to the nature of the heat source which 

in actuality has a very definite shape. Nevertheless, 
Rosenthal's initial attempt provided the impetus which spanned 


a succession of attempts to include all the important varia- 


bles in one well organized model. 


Grosh, Hawkins, and Peavane. investigated the effect of 
considering the thermal properties as a function of tempera- 
ture and obtained slightly better results. BESS modified 
Rosenthal's model by including possible heat losses to the 
environment but concluded that, at least for convection, 


the effect is negligible in comparison to the heat conducced 
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in the metal. 


The most important attempts to derive a satisfactory 
model concerned those which examined the shape of the heat 
Seunce. In 1952, tenes & proposed an equation from which 
the length and width of the weld pool could be determined 
as a function of the heat input. Subsequent experiments 
indicated that the relationship was reasonably accurate. 
Christensen, Davies, and Gjermundsen* examined the shape 
of the weld pool and its relation to a wide variety of 
welding conditions and associated variables. While their 
objective was not to formulate a model, they did present 
experimental data from which viable conclusions were drawn: 
a) pool width and length increase with heat input for a 
given weld speed and plate thickness, b) the pool is wider 
than predicted by point source theory, and c) the pool has 
a distinct "tear-drop" shape. These basic observations per- 
Mitted the derivation of a very successful model of tempera- 
ture distribution during welding. 


13 and Ta inealente bias constructed and 


In 1967, Pavelic 
verified a model which combined all the disparate elements 
necessary to predict temperature variations in the haz. 
Unlike Rosenthal, they considered the Variatiom or conduc] 
tivity and diffusivity with temperature, a finite heat 


source, heat losses to the environment through convection 


and radiation, non-linearites in the heat flow equation, 
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and preheating of the metal from arc radiation. Their ap- 
proach to the derivation of the model for a tungsten-inert- 


Gas weldamg process formea the basis for this thesis. 


C. BOILING HEAT TRANSFER 


The literature on boiling heat transfer was investiga- 
ted selectively with respect to possible mcdels relating 
to sub-cooled pool boiling on horizontal plates at a non- 


uniform temperature. 


The complexity of the boiling process has prevented 
researchers from developing a single formulation covering 
the wide range of temperatures over which boiling will oc- 
cur. Figure 1-1, reported by Rohsenow, nS is considered to 
represent the regimes of boiling, each regime having a 


uniquely formulated model descriptive of heat transfer. 


Rohsenow~> has correlated the data of Addoms for nu- 
cleate boiling on platinum wire in water at saturation tem- 
perature. Although the data was given for a specific ge- 
ometry, Rohsenow observed that the same results applied to 


horizontal surfaces. 


Correlations in the transition regime have not been 
firmly developed. Jordan,*> reporting the data of Drew 
and Mueller, Farber and Scorah, and Braunleh, concluded that 
if the data were plotted on a log-log graph, they would lie 


along a straight line from the critical flux (or, burnout) 
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to the minimum film boiling heat flux. The slope of the 


line varied with the surface finish: the steepest slope oc- 


CUrnGwAga@ta th a marnor smeothefinish. 


Chang, Botushansky, and Berenson*> developed correla- 
tlememeoco film boiling heat flux from horizontal surfaces. 
Jordan*> has reported that the Borishansky relation fits 


the data for pool boiling in water. 


D. SUMMARY OF PROCEDURE 


The development of a model for predicting temperature 
changes in the heat affected zone of an underwater weld 
required the construction of a theoretical model followed 
by experimental verification. The model was to encompass 
the latest techniques, and should be reasonably accurate 
in determining peak temperatures and cooling rates. 


The derivation of the model followed the approach taken 


20 in their work for surface GTA 


by Pavelic?? and Tanbachuchi 
(Gas Tungsten Arc) welding. The two-dimensional equation of 
heat flow with external sources was reduced to quasi-steady 
state form and transformed to its finite difference equiva~ 
lent. The external heat sources were examined and the proper 
relations chosen. Then, with the help of Pavelic's results, 
the equation describing the shape of the molten weld pool 


was derived. The resulting derivation included numerical 


constants which were determined through experimentation. 





18 
Finally, the computer Pacman was written. The method for 
lavyingeeuememeomranite difference grid was that ee iorea by 
Tanbachuchi. His solution was seFimed by having the shape 
of pool internally generated instead of being supplied as 


input. 


The experiment was conducted in two parts. First, the 
correlation constants in the pool configuration equation 
were determined by measuring the dimensions of molten weld 
pools while varying the welding parameters: Heat input 
and welding speed. Then, temperatures occurring during 
welding were determined by laying weld beads on submerged 
plates in which thermocouples had been installed. The re- 
sults of the final experiment were compared with the values 


predicted by the model. 
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CHAPTER II 


MATHEMATICAL MODEL AND COMPUTATIONAL TECHNIQUE 
A. PRELIMINARY OBSERVATIONS 


Welding is a process of joining two or more materials 
together by applying sufficient quantities of heat to a lo- 
cal area at their interface such that melting occurs and a 
union is achieved. The amount of heat required is on the 
order of wae BTU/HR and, for structural applications, the 
manner of generation is by electric arc. The method for 
directing the arc is by an electrode which may (consumable) 
or may not (non-consumable) contribute metal to the finished 
bond. The electrode concentrates the arc and creates a small 
pool of molten metal which is typically "tadpole" or "rain- 
drop" shaped when viewed in the plane perpendicular to the 
arc axis. As the electrode is drawn along the material a 


bead of weld (i.e., a chain of overlapping pools) is formed 


which represents the physical evidence of the interface bond. 


The largest percentage of the heat developed by the arc 
goes into raising the temperature of the material to its 
melting point. The remainder is accounted for by conduction, 
convection, and radiation, joule heating, and pre-heating of 
the material in the vicinity of the arc. Pre-heating will 
occur in all but submerged arc processes and is due to the 
divergence of the arc which takes the appearance of a cone 


with the electrode at its apex. The central part of the 
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arc, or the plasma column, is responsible for the creation 
of the molten pool. The remainder of the arc, while not 
generating sufficient heat to melt the material, will cause 
presacaeamae Hence, heat may flow into the material direc- 


Piveeeemrne pool and from the arc by pre-heating. 


While the above heat balance has been observed in nor- 
mal sruface welding, its complete relevance to the under- 
water welding process is less certain. Conduction and radia- 
tion have importance, but convective losses must be modified 
by boiling heat transfer phenomena. Joule heating losses 
are usually estimated to account for 1-2% of machine power 
and, hence, may be neglected. A search of the literature 
did not reveal any information on the dynamics of an arc 
surrounded by water so that the efeect of pre-heating could 
not be determined. For the experiment, it was decided that 
the effect, if it did exist, could be reduced to negligibi- 
lity by the proper choice of electrodes. Consequently, an 
electrode that could be placed directly on the metal sur- 
face without shorting out was selected. These "drag" type 
electrodes operated in a manner such that the metal core 
melted faster than the flux covering the flux coat, which 
came in direct contact with the metal, formed a hollow bar- 
rel and prevented arc diffusion. Consequently, it was assumed 
that pre-heating was negligible and all heat flowed into 


the metal through the molten pool. Interaction with the 
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environment was assumed to occur by radiation and boiling 


heat transfer. 


An approximation to the real mele process was then 
seen to be a boundary-value problem in two-dimensions with 
negative source terms. The interior boundary in the field 
was viewed to be formed by the contour of the molten pool; 
the exterior boundaries were considered as the four edges 
of a rectangular surface. The pool contour temperature was 
known to be the melting temperature of the material with 
the edges at ambient. Radiation and boiling heat transfer 


were included as the negative sources. 


B. EQUATION OF HEAT FLOW 


The governing differential equations are presented in 
this section. In the sections to follow the equations are 
reduced to their finite difference equivalents, the boun- 
daries are defined, and the relations for the source terms 


are determined. 


The general equation of heat flow with source terms 
2S) 


pcDT 


Ve (K T) + W, = TG 


where k,?0, and c are the conductivity, density and specific 
heat of the material and may be a function of distance, time, 
and temperature. W; represents the source terms and has the 
units, BTU/PT-HR. In partial differential notation for two 


dimensions, the equation becomes, 
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The equation in this form is for a coordinate system 
fixed to the control surface defining the temperature field. 
Since the interior boundary (i.e., the pool contour) moves 
through the field with velocity, V, the velocity of the elec- 
trode, the temperature at a given point will change with 
time. To eliminate the time dependence, the coordinate sys- 
tem was transferred to the moving boundary with the follow- 


ing transformation equations: 


Kk = x-Vt 

€ = ¢t 
where £ and € arc coordinates in the new system. (The de- 
tails of the calculation appear in Appendix 1.) With the 


transformation, the equation of heat flow became, 


The temperatures with respect to the new system could 
now be considered constant and were viewed to form isother- 
mal contours in the plane of the field with the pool contour 
forming the innermost. Hence, the pool boundary could be 


~“ 


called the isotherm of melting. 


The equation was put in non-dimensional form by 
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defining normalizing relations. It was decided to make the 
exterior boundary conditions homogeneous and the conditions 
along the pool contour unity. Hence, the normalized tempera- 


ture was defined, 


where Ty is the melting temperature of the material (ap- 
proximately 2700°F for mild steel) and sta the ambient 


temperature. 
The normalized dimensions were defined to be, 


as x ae Y 
x* = TaN * = a e 
Xx : Y xX 
oO Oo 


when these definitions were substituted back into equation 


2.1, it was found that Xo should be defined as 


x, = aj /V ‘ 


where @ refers to a specific value of the thermal diffu- 
Sivity. From a stability analysis of equation 2.1, recom- 


mended by Ames? 1t was determined that 


wiere aMin yvepresents the minimum thermal diffusivity for 


\ 


the material. Finally, the normalized dimensions became 
x* = 5 aa y* = —— Pie 
Min OMin 
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by cUrerzeUeton Of cquatlOns 2.2 and 2.3 into equation 2.1, 


the heat flow relation was obtained in non-dimensional form, 


3) ae 1 dk 4 oMin? zit. oT 
Jx7 7 oy- * k at a 3x’ +t al : re -T “= ox 

Fat 
where the "Star" superscript has been dropped. (See Appendix 


4 for details of the derivation.) 


Be aN th Dib PRRENCE SOLUTION 


Equation 2.4 was then converted to its finite difference 
equivalent so as to make it amenable to computer solution. 
Following the procedure recommended by Greenspar? for el dite 
non-linear partial differential equations. A finite differ- 
ence grid node (Figure 2-9) was defined and the bounding tem- 


peratures were approximated by a Taylor expansion about T 


Z 2 
T, =T, ©) ee at io ae i =1,2,3,4. 
OX 2 ax? 


Now, in the limit as the S; (i.e., the grid spacing) approach 
zero the variation in temperature between the bounding tem- 
peratures along a given axis will represent the actual rate 
of change of temperature at the node. For unegual grid spac- 


ing, these changes were defined as: 


aT 
On 
je ys Soeijee e SiS" Ho%5 
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The weighing coefficients were found by comparison with 
Taylor's expansion through combination of like terms. The 


second partials were similarly defined 


and the coefficients found from Taylor's expansion. Substi- 
tution of equations 2.5 and 2.6 into 2.4 resulted in the finite 


difference equivalent for the equation of heat flow: 


2a + S 


Tie \T) Dig pies BO ely 2 
sS-1S-FE— _ Sel 8 SS Sy 
as) S,+S,) aS,(S,+S,) S3 S,+S,) 


2 
4 e 
T 2 OmMin 
+ 4 | sce |* vie ym com \ [W. ] 
S,(S,+S,) KV (Ty T,) 





Sa S S 
iL Gls 12 2 T i a 
+ = an|| [¢ ) To + (eecgszy— ot erqsrtszy? | 
x aT S458, O S,(S,+S,) S5(S,+S5) 





/B 


2 
S55 S S 
3 Saat 4 3 
7 PONE a ey ere ome ee 
Shee S,(S3+S,) 3 S,(S3+S,) a | 
Sf 


2 


2 2 S978) 
a(Sjs 


) 
9) 





The details of the derivation are presented in Appendix 
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EoWatemeomez./ Operated on a grid (i.e., an assemblage of 

nodes) placed over the temperature field by assigning to each 
node point a value consisting of the weighted "average" of 
themes DOuUnding temperatures. Since the exterior points 
in the field were at a normalized value of zero it was de- 
cided that calculations should begin at the pool contour 
(T - 1) and move outward. A series of iterations over the 


grid were performed until the desired accuracy was obtained. 


C. THE BOUNDARY CONDITIONS 


This section deals with the derivation of the equation 
describing the shape of the pool contour (i.e., the interior 
boundary in the temperature field) and a determination of the 
location of the exterior boundaries. The treatment of the 
problem follows the approach taken by Pavelic?? in his study 
of the tungsten inert gas process for normal air welding. 
Pavelic's method for finding the shape of the pool can be 


summarized as follows: 


1. Rosenthal's two-dimensional equation with a mouine 
source of heat 1S analyzed to obtain the relations 
predicting the maximum dimensions of the pool: 
length and width. 

2. These theoretical relations are correlated with 


experimentally-obtained data. 
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3. The contours of actual pools are measured and a 
curve fitted equation is obtained. This equation 


is a function of the maximum pool dimensions. 


Hence, for any given set of welding variables, Pavelic could 
predict the contour of the pool. The details of Pavelic's 


analysis have been summarized in Appendix 3. 


It was found that the maximum theoretical dimensions 


Of pool are 


K @ ) 
; _ 7 Oem 
Ymax ~ ™m 1 ea 


For the maximum half width, and 


for the maximum length. The argument of the Bessel functions, 
Ky and Ky, is the radial distance from the origin of the arbi- 
trary coordinate system to the pool contour and was found to 


be a function of the lumped parameter, p: 


2 Seema 
27 a L A 
The relationship of Xmnax! Moar, and i to the pool is pre- 


sented in Figure 2.1. When plotted on a log-log graph, the 


pool width (= 2y ) vs p and x VS p appear as a straight 


max ax 


line and, hence, experimental data could be correlated by 
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the relations: 


2 =a » 
Ymax  ?P 


x = 2 
max oe 


with a, b, c, and d becoming the correlation constants. 


Experiments were conducted therefore to find values 
for the constants. Pools were made by welding beads on 
thin rectangular plates. At a predetermined point in the 
welding process the molten metal was ejected from the pool 
by a jet of high pressure air leaving a crater whose edges 
formed the desired contour. Only those craters nearest 
the ideal two-dimensional configuration were retained for 
measurement. An ideal 2-D crater has been defined in 
Figure 2.2. The maximum length and width were measured 
with a Bausch and Lomb microscope containing a movable 
micrometer table. A computer program was then written to 
fit a curve through the data points from which values for 
the correlation constants were obtained. The program with 


sample input is presented in Appendix 7. 


Two major problems were recognized early in the attempt 
to produce craters for measurement: a) a search of the litera- 
ture did not reveal any similar experiments conducted with 
stick electrodes, and b) it was very difficult to time the 
release of air to eject the metal during the under water 
welding process. Without the certainty that the metal was 


being ejected at the proper instant it was considered 
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possible that the craters might be distorted and measure- 
ments valueless. The lack of prior experiments conducted 
in air denied a basis for comparison to determine if the 
underwater data were reasonable. Consequently, it was de- 
- cided to conduct an initial series of runs in air where 
crater production could be closely controlled, merely to 


serve as a crude reference for the underwater results. 


Plates of two different thicknesses, 1/8" and 1/4", 
were used and the data plotted as shown in Figures 2.3 and 
2.4. The results for the 1/8" plate appeared to be in 
good agreement with the predicted relation, but those for 
the 1/4" plates were not. The lack of alignment was attri- 
buted to the failure to obtain two-dimensional pools in 
the larger specimens. The available power was insufficient 
tO Greate the required configuration, and, as a result, the 
crater formed was only a hollow depression on the surface of 
the specimen. Significantly, an increase in power (ExI) 
would have resulted ina larger value for the lumped parame- 
ter, p, which would have shifted the data points to the 
right and brought them more in line with the predicted 
relationship. (It was noted, however, that a line drawn 
through the 1/4" data points was approximately parallel to 
the theoretical line. This relationship suggested that the 
penetration of the arc could also be correlated with the 
lumped parameter, p, so that the equation for a 3-D pool 
INight eventually be derived. No attempt was made to do so, 


however.) From the plot.for air welding it was observed 
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that a) 2-D pools would not be obtained in water for 1/4" 
thick specimens and b) the width of the underwater pool 
should be larger than predicted while the length should 


be smaller. 


Pools were then created under water, measurements taken 
and the data plotted as shown in Figures 2.5 and 2.6. The 
location of the data points with respect to the theoretical 
line was noted and compared with air welded results. Since 
the location was as expected it was assumed that the measure- 
ments were reasonable. Curves were fitted through the data 
points. These curves were not found to be parallel either 
to the predicted or the air welded Laneee and the reason for 
the divergence (considered more fully in Chapter IV) could 
not be adequately explained due to the small size of the 
data sample and data scatter. Nevertheless, the correlation 
constants obtained from the fitted curve were used. It was 
consequently assumed that the length and width of an under- 
water pool were related to the lumped parameter, p, ac- 


cording to the equations: 


ihe (alt 
kos 
xu 4.42 p 
0.67 
ko = 
Wa SoG p 
where Sal and We are dimensionless quantities and 
a 24x% 
M M~Min 2 nies 
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with the maximum pool dimension relations defined it was now 
possible to derive an equation that would predict the com- 
plete contour of the pool. Pavelic suggested curve fitting 


an equation of the form: 


¥* = Bx*(x_ - x)? 2.8 


to which the boundary conditions were applied, 


(ly y =o @ 2a 

(2) yH=o @ a 

(3) Ye=o @ y= Wi’ 2 , X= KX, 
(4) Yy=o @ y=o,x=o0o0 


where x is the distance to the point of maximum width. 
After the applying of the boundary conditions, it was found 


that 


(W_/2)* (a+b) 27° 


eal (ax )@ (bx_)® 


where a and b are the correlating constants. The details of 


the derivation have been included in Appendix 3. 


Values for the constants were obtained by measuring 
the actual pool contours and fitting the resulting measure- 
ments with equation 2.8. A computer program, employing an 
M.I.T. non-linear least squares curve fitting routine, was 
written and has been presented in Appendix 7. It was hoped 
that a large number of pools would be available for measure- 
ment in order that the calculated correlating constants 


would be generally applicable. While a large sample size 
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was obtained for air welded pools, only one under water pool 
of acceptable quality was produced. Due to the accelerated 
solidification of the molten metal in water and crater walls 
took on a jagged appearance which made measurements diffi- 
cult although the maximum dimensions were generally easily 
determined. Consequently, it was decided to base the cal- 
culation of the constants on data taken from the air welded 
pools and then compare the predicted contour with the mea- 
Surements from the single under water pool. The result of 
this approach is presented in Figure 2.7. The predicted 
contour gave a wider head and narrower tail than was mea- 
sured but the overall shape was considered satisfactory. 
Nevertheless, it was realized that the contour equation was 
based on a sample size of one (i.e., the single under water 
pool). However, it seemed reasonable to assume that the 
general shape of a pool would be independent of the atmos- 
phere in which it was created site the dimensions were most 
strongly influenced by heat input and welding speed. With 


this assumption, the correlation constants 


were used to predict the shape of the under water pool con- 


tours. 
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De ~DERIVALTEONT Or THE HEAT LOSS TERMS, We 


For this model heat losses were assumed to occur 
through boiling and radiation. Although convective heat 
flux could be considered, it was estimated that these losses 
‘would account for only 3% of the total heat input to the 
specimen (15,000 to 20,000 BTU/Hr) and were therefore neg- 
lected. Consequently, heat losses were considered accord- 


ing to three distinct regimes as depicted in Figure 2.8. 


Nucleate boiling is characterized by the termination 
of bubbles at the heated surface which grow in size with 
increasing temperature. The agitation of the liquid near 
the surface caused by rubble activity -- formation, collapse, 
or ascent -~ Maintains a circulation of relatively cooler 
liquid near the surface which can result in heat fluxes on 


6 2 


the order of 10° BTU/FT“-HR. 


When the surface temperature reaches a value 50-90°F 

- above saturation the number of bubbles forming has increased 
to the point where they begin to overlap and create a solid 
film of vapor. The vapor film is not immediately stable and 
may collapse back to the heated surface where the process 
repeats itself. In this transition region the cooler liguid 
is brought into contact with the surface only intermittently 
and heat flux rates are not as high as they would be during 
nucleation. As temperature increases the film grows more 
stable, insulating the surface and reducing heat transfer. 


Finally, the film becomes fully stable and the heat flux 
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reaches its minimum value. At higher temperatures heat loss 


through radiation becomes the predominant term. 


The complexity of the entire process required that each 
regime be treated separately in determining the appropriate 
relation for calculating heat flux. First, however, it was 
necessary to define the temperatures which bounded each re- 
gime. Since nucleation begins near saturation temperature, 
the lower bound was chosen to be 212°F. Burnout, or the end 
of nucleation, will occur at a temperature which depends 
upon the condition of the heated surface, the amount of dis- 
solved gases in the liquid, and the ambient pressure. If 
the solubility of air at 14.7 psia is 6m 28/2, then, from 
data reported by Rohsenow™> burnout will occur neat 260°F. 
However, a mildly oxidized surface increases the tempera- 
ture, by an unknown amount. Consequently, 300°F (as the 
nearest round number) was defined as the temperature at 
which nucleation terminates and the transition regime begins. 
For the minimum heat flux, the data of Borishansky and Beren- 
son as reported by Joréan*> indicated that transition ends 


near 500°F. 


Once the boundary temperatures for each boiling regime 
were established it was necessary to determine the suitable 
heat flux relations. Rohsenow's*? correlation cf Addoms' 
data suggested that heat flux in the nucleate boiling regime 
varied according to the cube of the temperature difference, 


Ty 7 ‘sar’ 





a) 
The properties of the liquid and vapor were calculated at 


saturation and the proper relation was found to be 


ee eta en 


qn = [1.542 (Tb -Toam 


The details of the calculation have been presented in Ap- 
pendix 2. As indicated by deneane?, the data for film boil- 
ing on a horizontal surface in water is predicted withrea- 
sonable accuracy by Borishansky's relation. In this regime 


the heat flux was found to be 


2 


> a BTU/hr-ft ; 


dp = 32.4(T - Loam) ; 


where calculations have been summarized in Appendix 2. 


Heat flux during transition boiling has not been well 
defined. However, data reported in the literature indicates 
that measures fluxes tend to fall along a straight line 
plotted with log-log coordinates. This relationship sug- 
gested the derivation of an equation with the burnout and 
minimum film boiling heat fluxes as boundary conditions. As 
summarized in Appendix 2, the results of the derivation were 


found to be 


282 a 
) 


q. = 2265.11 silat 
z (Da sar 


where 


300° = Te = 500°F 





yy 
The equation for radiation heat transfer was determined 
by assuming the surface of specimen exchanges heat with the 
surrounding water. Hence, 
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Gp = 0.1317 x 10° , BTU/hr-ft 


B. COMPULADIONAL TECHNIQUE 


The computer program written to calculate the temperature 
distribution consisted of the four main operations: 
1. The shape of the pool contour is generated. 
2. The grid is layed out to form a computational block. 
3. The grid nodes are initialized using Rosenthal's 
temperature distribution equation. 
4. Successive iterations are carried out over the 


grid until the desired accuracy is obtained. 


A listing of the program is presented in Appendix 7, and 


instructions on its use in Appendix 6. 


Pool Contour Generation 


A value for the lumped parameter, p, is calculated from 
the input data and the maximum dimensions of the pool are 
determined according to the correlation derived in Chapter II. 
Then, using the contour equation, which is a function of the 
maximum dimensions, the computer generates a series of points 
which define the shape of the interior boundary in the tem- 
perature field. If projected on the x-axis (i.e., the 


longitudinal center line) of the pool, the points would be 
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Seaeede>/ 10080 an aGneh apart. This spacing Mallows from 10 
(0.5 inch pool length) to 20 points in defining the shape of 


the boundary. 
Grid Layout 


Except for minor alterations, the method of arranging 
the grid is identical to approach developed by Tanbakuchi~? 
and is sized approximately 100 lines to the inch in a recta- 
linear array. The grid is configured such that the specified 
thermocouple location appears at the intersection of a row 
and a column. A schematic of the grid forming a computa- 
tional block is presented in Figure 2.9. Columns are ar- 
ranged parallel to the center line (1.e., column number 1) 

of the pool; rows perpendicular. The intersection of column 
l and row 1 mark the origin of the coordinate system. In 


the output, rows forward, or in the direction of welding, 


are identified as positive quantities. 
Initializing the Grid Nodes 


Nodes within and on the pool boundary are given values 
of one (T* = 1) while the exterior boundary nodes are set to 
zero. As suggested by Tanbakuchi*” the intervening nodes are 
initialized with Rosenthal's temperature distribution equa- 
tion. Initializing the nodes considerably reduces the number 


of iterations required to define the temperature field. 
Calculation of Temperatures 


A Gauss-Seidel type iteration, which makes use of the 
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4? 
latest available temperatures, is performed until tthe speci- 
fied accuracy is reached. The program begins each iteration 
cycle at the interior boundary to avoid excessive computation 
in that part of the field initially at zero normalized- 
temperature. Calculation of the value for the heat loss term 
is accomplished by a search routine which locates the node 
temperature on the boiling curve and computes the appropriate 


Guame icy 
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CHAPTER III 


UNDER WATER EXPERIMENT 


The experimental verification portion of the thesis 
consisted of three steps: Production and measurement of 
molten pools made in air, production and measurement of 
molten pools made under water, and measurement of tempera- 


ture changes in a specimen being welded under water. 


A. EQUIPMENT 


The same equipment shown schematically in Fig. 3-1, 
was used for all three steps mentioned above and consisted 
of the following: 

1. Water Tank 

2. Gravity Welding Stand 

3. Power Source 

4. Blow-Out System 

5. Temperature Recorders 


6. Current, Voltage, and Time Recorders. 


Water Tank. The welding experiments were conducted in a steel 
tank measuring four feet long by three feet high by two feet 
deep. The exterior of the tank was sheathed in wood andwawaes 
ber matting to prevent electrical shocks. A platform, con- 
taining clamps for securing the specimen, was fitted to the 
inside of the tank. The platform was located one foot from 
the top of the tank so that the electrode, but not the elec- 


trode holder, would be submerged in the water. The tank was 
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filled by a line supplying ordinary tap water (Boston, 
Mass.) and drained by gravity discharge. The tank was grounded 


to the power source. 


Gravity Welding Stand. A Japanese gravity welder was modi- 
fied and fitted to the top of the tank. The device consisted 
essentially of a common stick electrode holder fastened to a 
carriage which could be moved along an inclined rod. A guide 
was fitted to the rod so that the electrode would draw a 
straight bead on the specimen. The holder was adjusted to 
the carriage until the electrode maintained a constant angle 
of fifty-one degrees relative to the surface of the specimen 
in the plane of travel and ninety degrees in the plane per- 
pendicular to direction of electrode travel. The tripping 
mechanism on the carriage, which rotated the electrode up 

and away from the specimen at the end of travel, was adjusted 
so that the arc would be broken after approximately six in- 


ches of weld bead was drawn. 


Power Source. A Lincoln arc welder capable of supplying up 
to 40 volts and 350 amperes direct current was used as the 
power source. The device was set to operate on reverse po~ 
faeLoye hires aebangememe was SelLeCctcad as Oeil > limecemsem 


power sources for actual under water welding. 


Pool Blowout System. A system for ejecting the molten metal 
from the specimen to form a weld crater was constructed and 


consisted of a nozzle, foot switch, and nose connected to a 
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400 psi air line. The nozzle was mounted on a swivel-necked 
holder so that adjustments could be made in directing the 
flow of air into the molten pool. The holder and nozzle 
combination were clamped to the welding platform inside the 
tank. Several preliminary experiments were conducted to de- 
termine the direction and angle of the jet. After consider- 
ing several possible locations, it was decided that having 
the nozzle discharge along the axis of the bead but opposite 
to the direction of welding would give satisfactory results. 
The nozzle angle was adjusted so that the air would jet into 
the head of the pool, directing the molten metal back into 
the tank, away from the operator stationed behind the gravity 
welding stand. Once the proper relationships were determined, 
the nozzle was securely fastened in place for the duration 


of the experiment. 


Temperature Recorders. Three clevite brush, Model 220, re- 
cording voltmeters were used to measure the output signal 
from the thermocouples. Each instrument contained two chan- 
nels so that a total of six thermocouples could be recorded 
Simultaneously. The event markers were connected electri- 
cally and repressing the key on any instrument would cause 

a Simultaneous voltage signal to be recorded on all six chan- 
nels. The instruments were calibrated at 200°, 700°, and 
1000°F by inserting a thermocouple, of the same type and 
dimensions used in the experiment, into a controlled bath. 


The chart speed was adjusted to 25MM/sec and the gain to 
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IMV/Division (or, 50 MV total scale). 


Common mercury-in-glass thermometers were selected to 
determine water and air temperatures. One thermometer was 
taped near the terminals of the recorders to measure the 
cold junction temperature; the other was suspended in the 
water to the depth of the specimen (approximately six inches 


below the surface at the far end of the tank). 


Current, Voltage and Time Recorders. A General Electric 
D-C (0.400 amp.), Model 10H, recording ammeter was used to 
measure the current flowing in the welding circuit. An 
Esterline Angus (0.100V) recording voltmeter was selected 
to measure the operating voltage during welding. A timing 
device consisting of a General Time Co., Model 320-1, elec- 
tric clock connected to a circuit breaker was used to mea- 
sure the amount of time required to lay a given length of 
weld bead. When the arc was initiated the breaker would 
close and start the clock; breaking the arc would cause the 
clock to stop.Dividing the length of the bead by the elapsed 


time thus recorded served to indicate the welding speed. 


All three instruments had been calibrated by the in- 
strument lab at the Boston Naval Shipyard and were assumed 


to be recording properly. 


B. MATERIALS 


The criteria used in selecting materials for the 
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experiment was based upon their being readily available at 
low cost and having wide usage in the industry. Three main 
types of materials were used: Steel plate, welding rod, 


and thermocouple wire. 


Steel Plate. All specimens were cut from cold rolled mild 
steel (A.I.S.I. 1010) sheets. Three different thicknesses 
were selected: 1/16", 1/8", and 1/4". For the pool measure- 
ment portion of the experiments the sheets were cut into 4" 
by 6" sections; the specimens in which temperatures would be 
measured were sized 6" by 12". Both sides of all specimens 


were cleaned by grit blasting to provide a uniform surface. 


Welding Rod. E7018-Al, 3/16" Dia., stick electrodes were 
found to be suitable for the experiment. This metal core of 
rod tended to melt faster than the flux covering disintegrated 
permitting the electrode to be "dragged" across the specimen 
during the under water operation. This feature was critical 
Since the arc would not be observed while under water. In 
essence, the flux covering formed a hollow cylinder which 
shielded the arc from the water, and maintained the proper 
length necessary to sustain the arc. Since the duration of 
submergence in water was short (20 to 30 sec!), the electrodes 


were not coated with any water-proofing materials. 


Thermocouples. Chromel-Alumel Type K, .020" Dia., thermo- 


couples were selected to measure temperature changes occurring 
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during welding. The selection of the size required a com- 
promise between the ability to withstand high temperatures 
(large diameters), and fast response (small diameters). 
Furthermore, the thermocouple junction should form a point 
(ideally) in the specimen to avoid thermal gradients. 
Consequently, the diameter chosen was the smallest considered 


practical for the experiment. 


C. EXPERIMENTAL PROCEDURE 


The experiment consisted of three parts: a) production 
of pool craters, b) temperature measurement, and c) reduction 


of data. 


Production of Pool Craters. Molten pools were first produced 
by welding in air, then under water. The essential steps for 


this part of the experiment were the same for both processes: 


1. Specimen aligned on work table. 

la. If an under water pool was desired the tank was then 
flooded until the specimen was covered by six in- 
ches of water. 

2. The power source was turned on and given a five 
minute warm up period. ~ 

3. An electrode was selected and inserted into the 
holder’. 

4. The holder and carriage on the welding stand were 


manipulated so that the tip of electrode was slowly 


brought to the point where the arc would start. 
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5. The carriage was drawn manually down the incline 
of the welding stand, laying a bead on the surface 
of the specimen. | 

6. The blow-out system was activated by depressing 
the foot switch at the same time the carriage 
reached the tripping point. 

7. The voltage, current, time, ambient temperature, 
and bead length were recorded for the run. 

8. The dimensions of the pool crater - maximum width 
length, and contour - were measured under a Bausch 
and Lomb microscope equipped with a movable, micro- 
meter table. For the experiment, the ranges of the 


variables involved were: 


Gagne tre 100 - 300 amps 
Voltage Zee SUV ores 
Welding Speed 50) = eZ b trina 
Ambient Temp 70 fae oS ac aee) 
40 - 43°F (water) 


Current was varied by adjusting the rheostat on the power 

source in steps of 25 amps., while the voltage was allowed 
to seek its equilibrium value. A particular welding speed 
could not be chosen prior to a run since manual control was 
required. However, an acceptable range of speeds were pro- 


duced by consciously moving the carriage along the guide 


rod. A mechanical system for better control of welding 





56 
speeds was attempted but the problem of starting the arc 


defeated the project. 


The problems connected with producing two-dimensional 
weld pools by stick electrodes were numerous and it is esti- 
mated that only 20% of the runs resulted in acceptable data. 
The welding stand did not allow the freedom of movement 
that is normally enjoyed in "Striking" an arc. The low volt- 
ages required so small a gap that the electrode was often 
grounded before the are could start. In addition, a précise 
relationship between the heat input and welding speed neces- 
sary to produce two-dimensional pools was observed. Since 
the speed could not be predetermined, there was no assur- 
ance that an acceptable pool would result from a run. In 
practice, a number of runs at various speeds were conducted 


at a given heat input until the proper pool was obtained. 


Temperature Measurement. All runs were conducted under water 


and the essential steps involved were: 


1. Preparation of the specimen. 

2. Aligning the specimen with the electrode. 

3. Flooding the tank until the level of water was six 
inches above the surface of the specimen. 

4, Turning on the power source and recording equipment. 

5. Welding a bead on the surface of the specimen. 

6. Recording temperature, time, voltage, current 


and bead length for the run. 
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The preparation of the specimen involved drilling holes 
at specified locations and mounting the Reese buclee. 
Twelve .040" Dia. holes were drilled in each plate in a man- 
ner that assured that the axis of the hole was perpendicular 
to the surface of the plate. (See Fig. 3.2 for relative lo- 
cations.) The centers of the holes formed a locus of points 
five degrees from the center line of the plate. An addi- 
tional hole (#12) was drilled near the edge of the specimen 
to record boundary temperatures. (During the experiment, 
however, it was observed that placing a thermocouple at this 
location -- about six inches from the center line -- was un- 
necessary since the temperature variation one-half inch 
from the center line was negligible!). The distance from 
the center line to the center of the holes varied from 0.125" 
for thermocouple locations #1 and #9 to 0.440" for location 
#8; the center-to-center spacing along the center line of 
the specimen was approximately one-half inch. This arrange- 
ment permitted temperature measurement within the heat 
affected zone but allowed sufficient spacing for mounting 


thermocouples. 


The thermocouples required special preparation to pre- 
vent their shorting out in the water. As shown in Fig. 3-3, 
water-proofing reguired sealing all joints with expoxy 
cement and covering the leads with a plastic sleeve shrunk 
to form a tight fit. (In their commercial form, the thermo- 


couples had already been covered with asbestos and silicon 
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impregnated fiberglass.) A ceramic bead, containing separate 
holes for each thermocouple, was installed to prevent the 
leads from making a junction outside the plate. The exposed 
leads were cleaned with emery cloth and inserted through one 
of the holes in the specimen. A pair of pliers were then 
used to draw the leads through the hole until the ceramic 
bead was nearly flush with the bottom of the plate. The 
edges of the hole on the top surface of the specimen were 
given a light peening with a small prick punch, and the leads 
extending above the specimen were then cut off and ground 
flush with the surface. The purpose of the peening was to 
give the thermocouples a preferential junction at the sur- 
face of the specimen, and avoid, as much as possible, thermal 
gradients in the specimen resulting from pools that were not 
exactly two-dimensional. In this way, six thermocouples 

were installed in each specimen, leaving six open holes. 
Twelve specimens were prepared such that no specimen dupli- 


cated another with respect to thermocouple locations. 


All specimens, with one exception, were 1/8" thick: 
the 1/16" plates burned through too easily and the penetration 
reguired for two-dimensional pools was not obtained with 
the 1/4" plates. One 1/4" specimen was prepared, however, 
in order to determine the response of the model to a situa- 
tion where the shape of the heat source was not the ideal 


confligu@vatwon. 


A prepared specimen, containing six thermocouples, 
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was then placed on the table in the tank, aligned with the 


electrode, and clamped firmly in position. (It was first 
hoped that the electrode could be made to follow exactly the 
center line of the plate, thereby simplifying measurements 
taken to locate the thermocouples. Actual runs revealed, 
however, that although the bead would run straight it could 
not be made to follow any predetermined path with sufficient 
accuracy. Conseguently, the positions of the thermocouples 
with respect to the longitudinal axis of the weld bead were 
re-calculated following the run.) The tank was filled with 
water until the specimens were submerged to a depth of six 
inches. The thermocouples were connected to the recorders 
and the power source turned on. The electrode was guided 
into position in preparation for starting the arc. At the 
instant of arc initiation the laboratory assistant depressed 
the event key on a recorder which placed marks on all six 
output tapes simultaneously. The carriage containing the 
electrode was then either allowed to move freely or was 
drawn by hand with a slight steady pressure down the incline 
of the welding stand until the tripping point was reached. 
The ambient temperatures (air and water), run time, current, 


voltage, and bead length were recorded. 


D. DATA REDUCTION 


The maximum dimensions of the weld pool were correlated 


to the lumped parameter, p, by means of a least squares curve 
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PReetno Lourie aval lablesas san Option from M.I.T.'s computer 
library. The details of the program have been presented in 
Appendix 7. The correlating parameters for the equation 
determining the contour of the weld pool were found by a 
non-linear least squares curve fitting routine. (Also see 


Appendix 7.) 


The Y-displacement of a thermocouple was measured along 
a line running through the center of the thermocouple hole 
and perpendicular to the center line of the weld bead. The 
X-displacement was measured along the center line from the 
point of arc initiation to the intersection of the perpendi- 
cular. A computer program was written to reduce the thermo- 
couple output. The program scaled the distances on the out- 
put tape by the ratio of the chart to the welding speed, 
and ensured that the measured temperatures were in the same 
coordinate system as the model. The program, with sample 


input data, appears in Appendix 7. 
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CHAPTER IV 


RESULTS AND CONCLUSIONS 
a. LEST RESULTS 


The experimental portion of the thesis consisted of 
two parts: a) determining the correlation constants used 
to define the shape of the heat source, and b) measuring 
the temperature distribution in the specimen during weld- 
ing. The results of the pool contour study have already 
been presented and described in Chapter II; it need only 
be added here that the results were in good agreement with 


theoretical predictions. 


The results of comparing the temperature distribution 
as predicted by the model with that measured by experiment 
is represented in Figures 4-1 and 4-2. The coordinates 
have been presented in dimensionless form. The ordinates 
May be transformed to real temperatures by the relation 


T = 2677T* + 503 , °R 


The abcissa for Figure 4-1 can be converted to inches 
with the relation 


x 


M2740 .59me , ie 


and for Figure 4-2 


x x*7 41.02) ae 


where the starred terms are the dimensionless values. 


In the figures, Yn 1s the distance from the weld center 
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line to the point of maximum pool width; Ys (Figure 4-1) and X40 
(Figure 4-2) are the locations in the specimen at which the 


temperatures were measured. 


Figure 4-3 has been included so that a comparison might 
be made between the predicted temperature distribution for 
underwater and air welding: the curve labeled "without boil- 


ing" represents air welding. 


Figure 4-4 presents the predicted temperature distribu- 


GIOn at Various GiStaenceswrrom mace Center finc. 


Figure 4-5 represents a comparison between measured 
and calculated cooling curve slopes as a function of distance 
from the maximum pool width. The slope of a given curve was 
measured at the average temperature and was defined to be 
the change in dimensionless temperature with respect to the 
change in dimensionless distance. For example, the average 
measured temperature for Run 104 (Figure 4-1) was found to 


be 0.4 and the slope at that point was -0.025. 


Ee ANALYSTS OF DATA 


Prediction of peak temperatures within the haz was found 
to be reasonably accurate. The average error was 13% of the 
measured temperature with the largest variation occurring at 


locations nearest the molten pool. 


The cooling portions of the predicted and measured curves 


were not found to be in good agreement. As a quantitative 
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indication slopes were measured at the average temperatures 
of the distribution curves and plotted as in Figure 4-5. 
Near the molten pool the predicted slopes were in error by 
as much as 200% and converged with the measured values only 


at the edges of the haz. 


Consequently, an evaluation of the effect of boiling 
heat transfer was made by neglecting such losses (i.e., as 
if the specimen were welded in air) with results presented 
in Figure 4-3. Such an analysis indicated that conventional 
boiling models were inappropriate when applied to the under- 
water welding process. However, nothing in the data suggested 


what mechanism might be controlling heat losses. 


It can be speculated, nevertheless, that the gases gen- 
erated from the oxidizing flux and possible dissociation of 
water molecules could disrupt the normal process of boiling 
heat transfer. These gases would collect at the tip of =n 
electrode to form a protective envelope over a portion of 
the specimen's top surface. Within the envelope the predomi- 
nant form of heat transfer would be radiation while exterior 
areas could support some boiling losses. The bottom surface 
(opposite to the side being welded) would not see the shield- 
ing gases and, consequently, heat losses would occur through 


boiling. This mixed form of heat losses would explain the 


slower rate of cooling observed during the experiment. 


Unfortunately, the process just described has not been 


described in the literature and the development of the 
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appropriate model was beyond the scope of this research 


study. 


C. CONCLUSIONS 


The object of this thesis was to develop a model that 
would predict the temperature distribution in thin plates 


welded underwater. 


It was found that the model was reasonably accurate 
in predicting peak temperatures within the heat-affected- 
zone, but that the degree of error in predicting cooling 
rates was unacceptably large. The initial assumption that 
heat losses can occur by conventional mechanisms of boiling 
heat transfer was, therefore, considered invalid. Rather, 
the proper model for heat losses must consider the effect 


of the gases generated during the welding operation. 
Given the above conclusions, it is recommended that 


1. Heat transfer fluxes be correlated for the 
case involving full boiling with radiation on 
one surface while the opposite surface supports 
Onl VeLadtariten and 

2. The study should consist of vertical as well as 


horizontal surfaces. 
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APPENDIX 1 


TRANSFORMATION TO QUASI-STEADY STATE 


Pn Egy, 


Ask aoa) keer dk oT eih en) 
fs +f5]/4+ Of a Se Ss — 
aE dy . OX OX ‘ dy oy " Ne Poze 


in the steady state system, the origin will be defined to 
be centered at the moving arc. Hence, the following trans-— 
formation-equations apply: 


x = X tv, 
t=¢€t 


where the primed quantities refer to the new system. 


Hence, 


0X = OX OX Y= Voce ey, ee ole 


xen 
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with these partial relationships, derivatives in the new 


system become, 
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Since temperature is constant in the steady state system 
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Substituting into equation (1) and dropping prime :-notation, 








ger o'r i fak ar , ak at], Mi _ -v ar 
ax2 ay2 k Ox OX dy oy k Q ox 
where 
a = k/pc 


The equation is simplified by noticing that 


1 ok oT __ 1 3k (aT) 
OX 


Finally, 


977 927 1 3k aT? oT 
Ree pee ie ES ea ers 
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APPENDIX 2 


PERTVATION OF HEAT LOSS TERMS 
A. Nucleate Boiling 


The Poneenows = COrrevatlLoOn als 


C.(T -T ) (g/A) gd 
iw USAT LL C _ UN ia ae eG be 7) 


he, Sfluphe, g (p71 -P,,) 


All quantities are taken at saturation (14.7 psia, 212°F) 


from the CRC HANDBOOK OF TABLES FOR APPLIED ENGINEERING 


SCIENCE. 
oF = 1.005 BTU/1b-°R 
Pry ae) ee 
ST te BAS ps 10°’ 1b-sec/£t? 
G9/g = 1 
he = 970.3 BTU/1b 
Py = Soe lib, fee 
oe = 0,087 Ieee,” 
9 = 40,4 % 10 1b/£t 


Substitution of these values into the nucleate boiling 
equation yields, 


1 BTU nee 


qn = [1.542 (T. -Toam) 


in normalized form, 


3 
qx = [1.542T* (Ty -Toag) } 


where, 
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Be Ean Beoadking 


The Borishansky relation reported by Jordan“? is 
3 > st 1 
re 0.028 | ee yt Pal ee ee 
F u,[V/g(p,-p.) 12 w SAT 
where, 
Oy ESE Wey ea 
px = 0.0374 Wsy/ize* 
B= ood se Ngee as Aee 
heg= 970.3 BTU/1b 
By Se aS & Moo! 1S eeBeye2e2 
i MS Oe MoT eae ee 


Substitution of these values into the film boiling 
equation yields, 


D.7> 


= ae ae? 
tees O84 (GY an , BTU/hr-ft 


SAT 


where the normalized form becomes, 


qt = 32.4{(T*(T, - ieee 


Ee Tsar 


C. Transition Regime 


The proposed relation for heat flux in this regime is 


= 7 b 
Gp = M(T-Toaq) 


to which the boundary conditions are applied: 


i Gn = Gy at T, = BOO. Th 
2) Gp = Ap at TO = 500°F 
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Baen, 


3 
[1.542 (300-212)] = 2.34 x 10 
0.75_ 


la) MReOo-212)° 


2a) M(500-212)? = 32.4(500-212) 2265.11 


for which 


m 2265. 1a 
b= 5.86 


Now, the transition heat flux equation becomes, 


By eis 


288 
) 


Dee 


Ap =re2 76951 
nd Ssy ANSE 


and the normalized relation, 


Bi oe 
288 
eal = 2265.11 Teer Ln 
: a (Ty Tsar? 


D. Radiation Heat Transfer 


For the heat flux during radiation, 


¥ 4 4 
Gp = 0 POPLCT —Toan) - 


Since the specimen is completely enclosed by the liquid 


environment, 


and 


it 


12 
=) 


1/(1/e,tl/e, - 1) 
the emissivity of steel cleaned by sandblasting is approxi- 
mately, 


Qe) = 0.8 


while the water is very nearly a black body and 


eo = 0.95 
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The heat flux equation becomes, 
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and in normalized form, 


ko = * = = 
ap OSE cy SEO, cE ae Tean)* gan] Toni 





78 
APPENDIX 3 


Derivation @f Pool Boundary Equations 


Pavelic's derivation of the pool boundary equation can be 


summarized as follows: 


1. Rosenthal's point source equation is analyzed 
to obtain the maximum theoretical pool dimensions. 

2. A correlation of measured data is obtained based 
on the predicted relations. 

3. A two-dimensional equation is fitted through the 


boundaries of experimentally obtained pools. 


1. The maximum dimensions of the pool. 


Rosenthal's equation for the temperature distribution 


from a point source, 





T- 7, = —% xk (ivr) exp{-Jvx] 
27 LK 


where 
Q = 3.415 EI , BTU/HR 
Ko (Avr), Bessel Function 
4, Reciprocal of the average diffusivity 


Vee Welding Speed (Ft/Hr) 


L, Plate Thickness (Ft). 


With the normalized relations, 








(es 


The temperature equation becomes 


3.4 See 1 
2TaL H 
m 


Te = K(r*) exp [-x*] 


Define the lumped welding parameter, 


ee So. 415 en 
27a L H 
m 
where, 
Hn = oe (Tmt) 


so that (with star notation dropped), 


T= PK, (r) exp [-x] 





The maximum dimensions are found since, 


| 


a O at maximum width Lee 


r= x and y = 0 at maximum length. 2 


From condition 1., the value for the maximum width is 


es 1 - [2 : 
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where Lm fonpaeraned at ff —]—=1 ana 
P. GExp | eye Ge 


Ko (ri? 


PaeeoOoL COntour Equation: 


An eguation is arbitrarily selected to match the boun- 


dary conditions for the pool, 
Y=o a dat <a aa) ee 


y= 0 , at Y=" Ypy, 


y= 0 , at x = O, 2 





™m 
Reena aite is then, 
fe gg Sl ce) 
From the boundary Sondcion y =Oaty-= eer 
B= Yn 


| 


m Len n 
aw =e 


where "m" and "n" are parameters to be determined experi-~ 


mentally. 
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APPENDIX 4 


NORMALIZING THE HEAT FLOW EQUATION 


The objective is to nondimensionalize all quantities. 


Consequently, the dimensionless temperature is defined to 


be, 
foe Uk 
Te = Z 
tn ~ To 
where, 
To = Ambient temperature 
Tn = Melting temperature of the metal 
and, 


Ace , at the weld pool 

T* = 0 , at the specimen edge 
the coordinate dimensions are transformed by a length scale, 
H, s® that 

x* = x/H 

Vie YA 
By making the appropriate substitutions, we convert the 


equation of heat flow 


et eo eee ok aT 2 aT 2 it Ss ee 
jx2 * Dy +k OOT 3x? 7 + kp Cepia) ee 


to non-dimensional form, 


- 2m d JU i H* 2 -VH OT 
oe ee eS St AS) t+ HR -apagl= 2S 








~ 


X y aT y @ dx 


with the star notation dropped. 


A value for H can be obtained by examining the 
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convergence criterion* for the related linear equation*, 


o2y oo’ Her Gg gq VH 3t _ 
se t apr + —¢(iR-7B) + — 5x = 0 


ay 


Ames” suggests that the solution will converge when 

0 < S < 2a/VH 
where the value, S, is the grid spacing for the numerical 
solution. We arbitrarily establish the upper bound on §S 
to be one so that 

H = 2a/v 
Since a small grid spacing is desirable for accuracy, we 
allow the diffusivity, a, to assume its least value. 
Hence, 


H = 20a /V 


Finally, the complete non-dimensional equation of heat flow 
becomes: 


977 9°T lel oT 
Im? * Oye ° k aT ce + { 





dy’ |* Lv2K ion nor 


2 
aT) tOin 2CMin 3T 
x y 
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APPENDIX 5 


FINITE DIFFERENCE EQUATION OF HEAT FLOW 


Consider the following grid node: 


1K (X > Y+S3) 
9 





4 
Tg(X, Y-Sq) 


The Ss. represent grid spacings between nodes and To 


will be a weighted "average" of the values at the four sur- 


rounding nodes. The following treatment is a summary of the 
approach employed by Tanbachuchi~? and presented in detail 
6 


by Greenspan. 


First, Taylor's expansion of the bounding temperatures 


about ps is taken: 








Si 977 
Ves opto tee kien (1) 
OX 2 1 oe 
oT S5 ACU: (2) 
T, = To7S2 ax ae > nee 
oT S3 9°77 (3) 
oT Si) 9°T (4) 
Ty = Don os + a 3y* + 
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where, 


a, = £ (Ss ) 


by comparing equation (5) with equation (1) and (2) and 


collecting like terms, we find 


 . S47S5 ee So - Sea 
aC = , Si Sac oe 
O Sy So M0 Sy S,+S,) “ So (Sy +S.) 
Similarly, for the second derivative at To 
o*T = at +a.,T, + awT (6) 
Ox ayer Abe. 2° 2 


where comparison with equation (1) and (2) yields 





a2 Z 2 
a= , a Faas ia ll Soe 
Oo S185 di S, (S,+S,) 2 S,(S,+S,) 


changes in the y-direction are treated in like manner. 
Substituting the different forms for 3T/dx,dT/oy, 

92T/9x7, andj?T/ gy2 in the equation of heat flow, we 

find the value of To as a function of Thy To, T3, Tyr and 


the grid spacing. 
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APPENDIX 6 


Notes on the use of the computer model 


Cost and Size 

The computer program is Fortran coded and consists of 
approximately 1300 input cards. Runs were done on an 
IBM 360 using the G compiler. The program requires be- 
tween 400 K and 450 K bytes of core and generally took 
six minutes to execute. With low priority specified, 


then, costs averaged around thirty dollars per run. 
Stability 


The model is data sensitive and will diverge rapidly 
with the slightest provocation. Stability can be con- 
trolled by input of the proper grid multipliers: 

Y MULT 1, Y MULT 2, Gx 12, and cx 34 . The model will 
not converge if these multipliers are greater than one. 
The particular values required for convergence depend 
on the magnitude of the welding speed: higher speeds 
require smaller multipliers. As a reference, it was 
found that a speed of 50 ft/hr and a multiplier of 

0.75 will assure convergence. Unfortunately, as the 
multipliers are decreased the core requirements in- 
crease to the point where the program exceeds capa- 
loblbshiesi22 


OU CDUE 


The output is presented in terms of node point temperatures. 
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Columns run parallel to the welding direction with 
column one identifying the long axis of the weld pool. 
Rows run perpendicular to the welding direction. The 
intersection of row one and column one mark the loca- 
tion of the origin of the system. Rows with positive 
incremental distances are read forward of row one in 
the direction of welding; these with negative distances 


in the opposite direction. 


The output may be plotted in two ways. Each node tem- 
perature may be transferred to a grid, thereby identify- 
ing the isothermal contours surrounding ere molten pool, 
or, the temperatures along a particular column may be so 
plotted to indicate how the temperature varies at some 
desired point from the center line of the pool (column 


one being at zero distance). 


Decemumcion OL Enput Cards 

1. First sixty cards are Bessel functions: argument, 
Ip: Ij, Ko, K,- FOrMat (on lcm 

2. NOPRP. Format (I3). Number of temperatures for 
the thermal properties of the specimen. 

8. TRED, TKRD, TDRD. Format” (32105) et etperatumer 
Conductivity, and Diffusivity Of eNom pee nen 


4, VOLT, AMP. Format (2F10.5). Veltage and amps tor 


particular welding run. 





0 . 


ales 


2% 


hee 


Lb ge 
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Diet KRD, DISTV, TIME. Format (4F10.5). Melt- 
ing Temp. OS ae Specimen Thickness (IN), are 
travel distance (In), Time (Min). 
TA. Format (F10.5). Ambient Temp. (°R). 
CRAD. Format (E12.4). Radiation Heat Transfer 


@eetiterent = o (E, +E BTU/f£t*-hr-R*. 


he 
AVCON. Format (F10.5). Average Conductivity of 
Specimen, BTU/ft-hr-R. 

XLONG, YLONG, XLONF. Format (3F10.5). Length (In) 
and Width (In) of Computational Block. XLONF is 

the distance in inches from front of pool to forward 
edge of block. 

Gvere, 6X345°  YMULT1L, YMUL¥2" Format (4F1l07o).% 

Grid multipliers which control the fineness of the 
grid: forward of the arc, behind the arc, from the 
base line to the thermocouple, and from the thermo- 
couple to the edge of block. 

NWANT. Format (I3). Specifies the number of grid 
lines between each two points on the molten pool 
boundary. (2 works well.) 

ITCHK. Format (I3). Number of iterations. 
Vee — 50): 

ITCHK. Format (I3). Number of iterations at which 
OUutDUE iS printed to check convergence. 


ICONV. Format (I3). Number of convergence checks. 





d Dis 


Joy 


17. 


18. 


dB JS 


Zr 


Zils 


Zee 


PRES 


24. 


25. 


26. 
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Pee Format (13). Number of iterations forward 
of the pool. 
NDONT. Format (13). Number of rows that need not be 
printed since convergence was not attained. 
PeERT. Format (£3). “1" if printed output 15 desuca, 
oO) Sie NOt . 
NEVIS1. Format (I3). First column to be printed. 
NEVIS2. Format (I3). Last column to be printed. 
NEVISI. Format (13). Steps between first and last 
column. 
EPSI. Format (F10.5). Tolerance for calculated 
Temperatures. (.001) 
AKHAR. Format (F10.5). Coefficient in Rosenthal's 
CqrariOn  fOn Initialezing grid. 
JWAN. Format (13). Number of curves to be plotted. 
(2 maximum) 
NTC. Format (13). Number of thermocouples or points, 
required (max of 2). 
Y¥TC(I). Format (F10.5). Distance in inches for each 
thermocouple from base line. 
ETCN(DT)> Formac €15)." Thermocouple number, ot any 


dummy identification. 
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APPENDIX 9 


EXPERIMENTAL DATA 


The experimental data has been arranged as follows: 

A. Data for the maximum dimensions of air welded 
pools. 

B. Data for the maximum dimensions of underwater 
pools. 

C. Pool contour measurements. 

D. Data for temperature distribution with sample 


thermocouple output tapes. 
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C. POOL CONTOUR MEASUREMENTS 


RUN #4 RUN #8 RUN #13 
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D. WELDING VARIABLES FOR TEMPERATURE DISTRIBUTION 


RUN# 
100* 
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THERMOCOUPLE OUTPUT TAPES 
CHART SPEED: 25 mm/sec 


VERTICAL SCALE: 1 MILLIVOLT/mm 
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